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We analyse the full counting statistics (FCS) of a superconducting junction in non-equilibrium 
in the limit of small interface transparency. In this Hmit we treat both supercurrent and multiple 
Andreev reflections on equal footing and show how to generalise previous results for both phenom- 
ena. Furthermore, we also allow for different gaps of both superconductors and investigate the 
intermediate regime which allows to make contact with previous results on normal-superconductor 
heterostructures. We also compare our predictions in this regime to experimental data. 
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INTRODUCTION 



The charge transfer statistics of superconducting tun- 
nel junctions have been an object of interest for many 
years.l^^ At energies below the gap two phenomena are 
observed: on the one hand, coherent supercurrent of 
electrons between the superconductors depending on the 
phase difference and on the other hand multiple Andreev 
reflection (MARs).'^^! Andreev reflection (AR) refers to 
the retroreflection of one electron from another super- 
conductor as a hole leaving a Cooper pair behind. Such 
processes may happen multiple times if the two leads at- 
tached are both superconductors. 

The properties of Andreev reflection and supercurrent in 
multiple hybrid systems have been investigated experi- 
mentally both for tunnel contact^l^ and more involved 
geometries involving quantum dots.'^ 
So far a complete analysis of the charge transfer statis- 
tics through superconducting junctions has been done 
either for the dc-componenlP or the ac-component 
in equilibrium.'^' Only the current has been anal ysed 
completelylSl but a complete discussion of nois^^^EIl and 
possible higher cumulants of the charge transport statis- 
tics so far is missing. 

In this work we want to calculate the full counting statis- 
tics (FCS), which allows to flnd the cumulant generating 
function (CGF) of charge transfer. This function allows 
to calculate all cumulants of the current flow. However, 
in order to be able to rationalize the result for the FCS we 
will restrict ourselves to the case of small transmission. 
Such restriction will allow for a reasonable truncation of 
the MAR processes as will be discussed in more detail 
below. This case will be presented in full detail, espe- 
cially also incorporating the possibility of different gaps 
for the two superconductors (SCs). 

In Section[ll]we will present the general model and frame- 
work we use to calculate the cumulant generating func- 
tion. Its actual form will be discussed in two steps: in 
Section III we will discuss the properties of supercurrent 



and in Section IV we will present the results for multiple 
Andreev reflections. For different gaps of the two super- 
conductors there is an intermediate voltage regime which 



will be analysed in Section |V| We will compare our re- 
sults to experimental data in Section [Vl] and conclude in 
Section Ivni 



II. FULL COUNTING STATISTICS 



We calculate the cumulant generatin g fun ction x(-^j 4') 
using the generalized Keldysh technique^^ES] which allows 
to proceed via the Hamiltonian formalism.'^ A (t) refers to 
the counting field and </> is the phase difference of the SCs. 
x(A, 0) in the case of charge transport is given by 



X 



(1) 



where Pq{4') is the probability for the charge q to be 
transferred through the system during a given (long) 
measurement time r. 

Partial derivatives of x(A, (f) with respect to A give direct 
access to the cumulants (irreducible moments). Mod- 
elling the contact between two SCs as illustrated in Fig. 
[ijis straightforwarcF' 

H{t) = Hl+Hr + Ht, (2) 
where {a = L, R) 

k,(T k 

are the Hamiltonians for the two SCs using units such 
that e = h = ks = I- The voltage is applied symmet- 
rically so that = —fJ-R = V/2 follows for the chem- 
ical potentials. The corresponding normal and anoma- 
lous Green's functions (GFs) are given in [14] and are 
abbreviated as gaai^) and /q(w), respectively. The 
quasiparticle DOS is strongly energy-dependent Pl/j^ = 

P0,L/fi|w|/^w2 _ /^2^^^ 

The tunnel Hamiltonian has to take into account the 
time-dependent phase of the SCs 



i+i?,+e-'*M/2^+i,), (3) 
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FIG. 1. Sketch of the system under investigation: two su- 
perconductors are tunnel coupled with a tunnel amplitude 
7. Additionally they have a time-dependent phase difference 
(pir) = (1)0 + 2eVr. 



where 7 is the amplitude of the tunneling coupling and 
Lo-, i?CT refer to the field components at a; = where 
tunneling is assumed to occur. 

To study the FCS we calculate the CGF as a general- 
ized Keldysh partition function. The connection to the 
system's Hamiltonian in Eq. H is given bjK^ 



lnx(A,0,T) = ( Tcexp 



(4) 



where T'^'*) denotes Eq. ([s]) after the substitution — 
i^o-e^*'*''*''^^- C means the Keldysh contour and Tq means 
time ordering on it. The measuring field X{t) has to be 
both time- and contour-dependent. It changes sign on 
the different branches of the Keldysh contour to account 
for the charge transfer. Additionally A(i) is nonzero only 
during the time interval [0,t]. We use the standard 
expressiorfi^ 



91nx(A,(^,r) 
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Tc 



dX 



(5) 



to find the CGF as the counting field derivative of T^. 
Compared to the case of tunnel contacts between normal 
metals and SC^^^Hi^ the counting field derivative in Eq. 
([5]) not only has normal but also anomalous contributions 
(due to superconductivity) leading to 



dX 
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2^ 



normal contribution 



anomalous contribution 



(6) 



where G and T refer to the exact-in-tunneling and A- 
dependent GFs and indeces ±± refer to the positions of 
the two time arguments on the Keldysh contour. 
Eq. ^ can then be integrated with respect to A to access 
the CGF. The normal contribution gives rise to MARs 



and quasiparticle tunneling, whereas the anomalous con- 
tribution gives rise to Josephson tunneling. We discuss 
both parts separately. 



III. JOSEPHSON TUNNELING 

We first evaluate the second part of the expression in 
Eq. ^, meaning 




duj 



2 2n 



As discussed before we can evaluate J-^ and J-ji exactly 
by means of their corresponding Dyson equation.!^ How- 
ever, the result is more transparent if we only give the 
first order in F = tt'^j'^ PolPor, keeping in mind that this 
is reasonable only for T <C ^l/r- In this case we can 
use the approximation of the GFs as in [50], where the 
tunneling self-energy is real and purely off-diagonal for 
energies below the gap and diagonal for energies above 
the gap. The result for the CGF is then 

lnXa,i(^>>T) 

— In { 1 + r„ [(e'^ (cos <?!> -I- sin 0) - 1) 
27r 

+ (e-*^(cos0-sin</>) - 1)]} 
' Al- . ( ^r-\ujl\ 



where 



A, 



(7) 



" (Ai-c.i)V2(A2j-^«)i/2' 
using ul/r =uj- hl/r- 

Of course, this result is in perfect accordance with the 
form in [1 , where, however, the transmission coefficient 
is slightly different as there the case of a multichannel 
diffusive contact is treated, whereas here we discuss the 
case of a single channel ballistic contact. 
Compared to Refs. |2II9| we only observe the first har- 
monic of the current in the Josephson frequency since 
we only treat the case of small interface transparency. 
Indeed, we immediately find for V < 2A]^/jj 

(/a,i)(T) cx sin[(f){T)], 

which corresponds to the dc- and ac- Josephson current 
depending on whether a voltage is applied or not. 
Calculating the noise from Eq. Q we find that it may 
become negative even when no voltage is applied. This 
peculiarity is well knowrP and can be traced back to the 
fact that probabilities Pq{4>) in Eq. ([T) can become neg- 
ative for SCs due to the additional aependence on the 
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phase difference 4>. In this case the interpretation of 
Pq{4>) from the Wigner representation as a reasonable 
probabiHty is rendered impossible. This fact is discussed 
in more detail in Refs. [2 [52]. 



IV. MULTIPLE ANDREEV REFLECTIONS 



We go over to the evaluation of the first part of Eq. 



dT, 



\{t)' 



dX 



T 
' 2 



E 



2^ 



(8) 



From the second part of Eq. ([6]) we obtained an ac- 
current contribution for V ^ 0. Eq. (|8| will now pro- 
duce further dc-current contributions known as MARs 
and discussed in |23I24| . Our calculation is closely re- 
lated to Ref. [S]. We first discuss V/2 < A^/^ and will 
complete our calculation later. 

The time-dependent coupling in Eq. ^ allows for a fi- 
nite result for Eq. ([s]) even for voltages below the gap. 
For the following calculation it is easier to introduce com- 
bined GFs of the normal and anomalous part 



L/R. 



We have used the abbreviations 



S,L/R = 7 



9n,L/R 



++ 

9n,L/R 9n,L/R 



fn,L/R. 

J 7 



n.L/R 



,-'iA f — h 



Jn,L/R 



L/R 



f++ 
Jn,L/R 



In the case of low transmission the electronic transport 
can be described by a sequential tunneling picture^ ^ 
with transmission coefficients given by a product of tun- 
neling rates for each MAR. The corresponding transmis- 
sion coefficients are given by 



a,n,L/ R 

'A 



4r" Yl 



^L/R 



fe=l ^L/R 



(w T keV/2)^ 



X0 



L/R. 



^TkeV/2\ \ LJHj'J 



A2 
^R/L 



^L/R , 

Aj^/i-|w±fceF/2|\ 



A 



R/L 



\uj-neV/2\\uj + neV/2\ 



^{uj + neV/2Y - /^\yj{u: - neVl2y - A\ 
f \u + neVl2\- Ar \ ^ ( \u - neV/2\ - Al 



Afl, J \ Al 

The result for the CGF for V < 2Al/r is then 



Due to the special time-dependence of the coupling ele- 
ments every GF admits a Fourier expansion of the form 



9L/R.it,t') 



^gm0(t')/2 



2^ 



which means that gL/R{^,^') = ^ + '^6V)5(e — 

e' -I- neV). Therefore, in order to calculate the different 
transport properties we have to find the Fourier compo- 
nents gn„i,L/R = 9nni,L/Ri^ + neV,uj+meV). The Dyson 
equation for the Fourier components is slightly more com- 
plicated as e.g. in the normal-superconductor cas^^^ due 
to the coupling of different Fourier components. We de- 
fine gn,L/R = 9lm{^ + neV), fn.L/R = fL,R{'^ + neV) 
and the Dyson equation can be expressed a^ 



G 



nm,L /B. 



9mn,L/R^n,m H~ ^nn^mrL.L/ 

+Vn,n-2G'^_2 m.,L/R ,n+2Gn+2m,L/Ri 



where 



Vn,n+2 
Vn,n-2 



^S,R/L9n~l,L/R ^s.R/L9n+l,L/R 
^S..R/L9n~l,L/R ^ s,R/ L9n+1,L/ R 

^SM/Lfn+1,L/R 
^SM/Lfn+l.R/L 

^S,R/Lfn-1,L/R 
^SJi/Lfn-i^L/R 



lnXa.2(A,T) 



duj 
2^ 



In l-f ^ Ta,n,a 
n=2, a=L/R 

X [(e'^" - l)nF(w - nV^/2)(l - nF{^ + nVl2)) 



-fe-*-^"nj.(w + nVl2){\ - npi^^ - nV/2))]} , (9) 
where 11^(0}) is the Fermi function. 

Of course, the results for the current and noise that can 
be obtained from Eq. ^ are identical to the ones ob- 
tained in [21 [TO] for tunnel junctions. 



V. ANDREEV REFLECTIONS FOR 
INTERMEDIATE VOLTAGES 

So far we have considered V < 2ALji in which case 
both SCs contribute to the MARs. However, for A^ ^ 
Aji there is an intermediate regime in which only one of 
the SCs contributes to the AR processes. The description 
of this intermediate regime is rather involved, however, it 
does not lead to surprising dependencies on the counting 
fields. We consider, without loss of generality, A^ > Afi 
and the following regimes for the energies 

• b: Ar <\lj + V/2\ < Al, \^ - V/2\ < Al/r, 

• c: Aii<\uj + V/2\ < Al, Ar < \lj - V/2\ < Al, 

• d: Ar<\lj + V/2\ < Al, \io - V/2\ > Al/r, 
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• e: \uj + V/2 < Al/r, Ar < \uj - V/2\ < Al, 

• f: |w + V/2\ > Al/r, Ar <\uj- V/2\ < Al, 

• g: \uj- V/2\ < Al/r, \uj + V/2\ > Al/r, 



lnXc(A,T) 



lnxd(A,T) 



lnXe(A,T) 



lnx/(A,T) 



lnxi(A,T) 



where we have used nL/R = np{uj ~ ^^L/R) and 
- (Ai - - Ar^u^I - A|)V2 - 

y ^ rl^Lll^fll ^ -p 

' 2K-Ai)i/2(c.|-A|)V2 - 

The other energy configurations do not contribute to 
charge transfer. We observe an interplay of AR and 
single-electron transmission with the typical dependen- 
cies on the counting fields. In this respect no new 
noise features apart from those already discussed in the 
literatur^i^! appear. 

The full result for the CGF for a superconducting tunnel 
junction have the form 

lnx(A,0,r) = lnXa,i(A>,r) -|-lnXa,2(A,T) +lnXc(A,T) 
+ lnXd(A,r) -|-lnXe(A,0, r) -|-lnx/(A,T) 
+ lnx»(A,T). (10) 



• h: + V/2\ < Al/r, \oj - V/2\ > A^/r, 

• i: Iw - V/2\ > Al/r, \io + V/2\ > Al/r. 

We considered each configuration and solved the Dyson 
series. Using Eq. ^ we obtained the following results 



I 

Of course, choosing Ar = we reproduce the result in 
Ref. |18| in the low transparency limit. 
Nonetheless, we want to show that the above derived 
expressions indeed allow for a correct description of a 
superconducting tunnel junction at finite but possibly 
different Al/r- 

VI. ZERO TEMPERATURE SIS JUNCTION 

We want to discuss how to reproduce the results in 
from the semiconductor model. We limit ourselves to 
T = and in the first part only consider A^ = Ar = A. 
Furthermore we limit ourselves to the lowest order in F 
regardless of the energy dependence so that the only con- 
tributing parts are lnxa,i and Inxi- This approximation 
is typically justified for the experimental situation of two 
superconductors coupled via an insulating barrier (SIS 
junction). We disregard the ac-component and are con- 
sequently left with the calculation of the current from 
In Xi ■ We arrive at the result also presented in ^6] 



= r y" ^ In {1 + AT,[nF{uJ ~ V){1 - npiuj + V)){e''^ - 1) + npiuj + V){1 - npiuj - V)){e-^'^ - 1)]} 
^2tJ ^\n{l + T,[nLil-nR){e'^-l)+nR{l-nL)ie-'^-l)]} 

/J 
^ In { 1 + Fa [(e^^ - 1) (cos <j) + sin </)) + (e^^^ - 1) (cos <j) - sin </>)] } 

^2tJ ^ln{l-H4F,K(l-nfl)(e'^-l)+nfl(l-ni)(e-^^-l)]} 
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( V-2A \ op rV-A doj ui{V-ui) T/ \ n 

V A / °^ JA 27r (w2_^2)[(y_^)2_^2]l/2 , V ^'J 

>|-2A \ op r\V\-A i^(\V\-Lo) V^n 

A i Ja 2Tr (w2-A2)[(|y|-a;)2-A2]i/2 ' ^ 



Since the two cases only differ in sign we choose V > from now on. Using elhptic integrals and reintroducing Sl-units 
we arrive at 

= (^^) { - ieV + 2A)[Kia) - E[a) 



where a = f^=|f and f = 'i-f'^n'^ PolPor.- Using K{0) = 
E(0) = tt/2 one obtains the well-known discontinuity in 
the tunneling current at eV — 2 A to 7r/4 times the tun- 
neling current in the case of normal contacts. 
The case of 7^ Aji cannot be treated integrated an- 
alytically anymore. We obtain for the current to lowest 
order in T 
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du 
2^ 



r|a;/j||a;L| 



2K-A|)V2(^2_A2)l/2 

X [uLil - ur) + niiil - ul)]} 
ujl\ - Al\ f\i^B\ - A_r' 



A, 



Ar 



(11) 



We compare this result to the experimental data from 
[?7] for a Nb-Sn junction in Fig. [2] 



A cusp is observed for the current at y = jA^ — A/jj and 
a steep rise is again observed at \V\ — Al -I- Ar. For 
\V\ = |Al — IS.r\ the DOS of electrons participating in 
the charge transfer in both SCs is maximal and allows 
for a calculation of the order parameters. 
We do not observe perfect agreement as the surfaces of 
the SCs in the SIS junction can be nonhomogeneous lead- 
ing to different characteristics of the SC at different posi- 
tions. Furthermore Nb is already a strongly coupled SC 
making BCS theory not perfectly applicable. Nonethe- 
less, the agreement is accetable and shows that the CGF 
in Eq. ( 10 1 allows for a description of the case A^, 7^ /S.r. 



3'^ 



'si 

(13 
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U 





FIG. 2. Current in a Nb-Sn junction at T = 3.38K with 



= 0.19 •10-'^ eV, A 



^ = 1.49-10"'^ eV. The experimen- 
is shown as the dashed curve and 



tal measurement from 
compared to the theoretical prediction using Eq 



(111, solid 



VII. CONCLUSIONS 



To conclude, we have derived the CGF for a supercon- 
ducting tunnel junction taking multiple Andreev reflec- 
tions, Josephson tunneling and possibly different gaps of 
the two superconductors into account. In this way our re- 
sult generalizes previous results for the CGF taking only 
certain aspects into account. The understanding of the 
CGF in simple tunnel junctions paves the way to under- 
standing more involved geometries. 
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esting discussions. 
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